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ABSTRACT 

We present a novel technique to determine the absolute inclination of single stars using multi- 
wavelength sub-milliarcsecond astrometry. The technique exploits the effect of gravity darkening, 
which causes a wavelength-dependent astrometric displacement parallel to a star's projected rotation 
axis. We find this effect is clearly detectable using SIM Lite for various giant stars and rapid rotators, 
and present detailed models for multiple systems using the reflux code. We also explore the multi- 
wavelength astrometric refiex motion induced by spots on single stars. We find that it should be 
possible to determine spot size, relative temperature, and some positional information for both giant 
and nearby main-sequence stars utilizing multi- wavelength SIM Lite data. This data will be extremely 
useful in stellar and exoplanet astrophysics, as well as supporting the primary SIM Lite mission through 
proper multi-wavelength calibration of the giant star astrometric reference frame, and reduction of 
noise introduced by starspots when searching for extrasolar planets. 
Subject headings: astrometry — stars: fundamental parameters 



1. INTRODUCTION 

SIM Lite is currently expected t o have ^80 spec- 
tral channels (|Davidson et al.l |2009[ ). spanning 450 to 
900 nm, thus allowing multi-wavelength microarcsec- 
ond astrometry, which no current or planned ground 
or space-based astrometric project, (GAIA, CHARA, 
VLT/PRIM A, etc.) is able to ma tch. We showed in our 
first paper (jCoughlin et al.|[2010f ). hereafter referred to 
as Paper I, the implications multi-wavelength microarc- 
second astrometry has for interacting binary systems. In 
this paper, we discuss an interesting effect we encoun- 
tered while modeling binary systems, namely that grav- 
ity darkening in stars produces a wavelength dependent 
astrometric offset from the center of mass that increases 
with decreasing wavelength. It is possible to use this ef- 
fect to derive both the inclination and gravity darkening 
exponent of a star in certain cases. 

Determining the absolute inclination of a given star 
has many practical applications. There is much interest 
in the formation of binary stars, where whether or not 
the spin axis of each star is aligned with the orbital axis 
provides in sight into the formation history of the system 
(jTurner et al. 1995). The mutual inclination between the 
stellar spin axes and orbital axis can greatly affect the 
rate of precession, which is used t o probe stel lar structure 
and test gene r al relativity (iSterne 1939a. bl3: lKopal|[l959l : 
lJeffervlll984) . lAlbrecht et al.l (|2009. ) recently reconciled 
a 30-year-old discrepancy between the observed and pre- 
dicted precession rate of DI Herculis through observa- 
tions which showed the stellar spin axes were nearly per- 
pendicular to the orbital axis. Along similar lines, extra- 
solar planets discovered via the radial velocity technique 
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only yield the plane tary mass as a function of the incli- 
nation of t he orbit ("Mavor fc Ouelojll995t iNoves et al.l 
1997; Mar cv fc But ler 2000), and thus, if one assumes 
the planetary orbit and stellar rotation axes are nearly 
parallel, determining the absolute inclination of the host 
star yields the absolute mass of the planet. If the stel- 
lar spin axis is found not to be parallel to the plane- 
tary orbital axis, this provides valuable insights into the 
plan et's formation, m igration, and tidal evolution histo- 
ries (|Winn et al.|[2006.: Fabry ck y fc Winn 2009 ). A final 
example is the study of whether or not the spin axes of 
stars in clusters are aligned, which both reveals insight 
into their formation processes, as well as significantly af- 
fects the determination of the distances to those clusters 
fjackson fc Jeffries 2010). 

Our proposed technique can also be used in conjunc- 
tion with other methods of determining stellar inclination 
to yield more precise incli nation values and other stel- 
lar parameters of i nterest. iGizon fc Solankil (|2003l ) and 
iBallot et al.] (|2006l ) have shown that one can derive the 
inclination of the rotation axis for a given star using the 
techniques of astroseismology given high-precision pho- 
tometry with continuous coverage over a long baseline, 
such as that provided by the CoRoT and Kepler mis- 
sions. This technique is sensitive to rotation rates as 
slow as th e Sun's, but b ecomes easier with faster rota- 
tion rates. iDomiciano d e Souza et al. (2004) discuss how 
spectro-interfcrometry can yield both the inclination an- 
gle and amount of differential rotation for a star, param- 
eterized by a. For both eclipsing binaries and transit- 
ing planets, the observation of the Rossiter-McLaughlin 
(RM) effect can yield the relative co-inchnation between 
the two components (iWinn et al.l l2006t lAlbrecht et al.l 
I2009t iFabrvckv fc Winnll2009D . The technique we pro- 
pose in this paper would be complementary to these tech- 
niques in several ways. First, it would provide an inde- 
pendent check on the derived inclination axis from each 
method, confirming or refuting the astroseismic models 
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and spectro-interferometric and RM techniques. Second, 
in principle the astroseismic technique is not dependent 
on the gravity darkening coefhcient /3i, and the spectro- 
interferometric technique is correlated with the value for 
a; combining techniques would yield direct and robust 
observationally determined values for i, a, and /3i. Fi- 
nally, the accurate, observational determination of a and 
/3i, (along with stellar limb-darkening), is critical to ac- 
curately deriving the co-inclination from the RM effect, 
as well as other quantities in stellar and exoplanet astro- 
physics. 

In this paper, we also present models for and dis- 
cuss the determination of spot location, temperature, 
and size on single stars, which produce a wavelength- 
dependent astrometric signature as they rotate in and 
out of view. Star spots are regions on the stellar surface 
where magnetic flux emerges from bipolar magnetic re- 
gions, which blocks convection and thus heat transport, 
effectively cooling the enclosed gas, and thus are fun- 
damental indicators of stellar magnetic activity and the 
internal dynamos that drive it. Il^ik et afl (|2007f ) discuss 
how the observation of spot location, duration, stability, 
and temperature can probe the stellar interior and con- 
strain models of magnetic flux transport. Through the 
observation of the rotation rates of starspots at vary- 
ing latitudes, one is abl e to derive the differe ntial ro- 
tation rate of the star (jCollier CameronI 120021) , which 
may be directly related to the frequency of starspot cy- 
cles. Mapping spots in binary star systems provides in- 
sight into the interaction between the magnetic fields 
of the t wo components, which can cause orbital period 
changes (IADDleg ate!ll992D , radii inflation (jLopez-MoralesI 
120071 : iMorales e t al. 2008), and may possibly explain the 
^^2-3 hour period gap in cataclysmic variable systems 
(jWatson ct al. 2007). Detecting and characterizing star 
spots via multi-wavelength astrometry would be com- 
plementary to other existi ng techniques, namely opti- 
cal interf erometry f Wittkow ski et al.ll2002f), tom ographic 
imaging ponati et al. 2006: Auriere et a l.ll2008D. photo - 
metric monitoring (Alcksccv 2004; Mo sser et al.ll2009f ). 
and in the future, microlcnsing (Hw ang fc Hanll2010( ). 

We present the details of our modeling code, reflux, 
in SJH discuss the inclination effect and present models for 
multiple stars in discuss the spot effects and present 
models in Sj4l and present our conclusions in ^ 

2. THE REFLUX CODE 

REFLU}(E] is a code that computes the fiux- weighted as- 
trometric reflex motion of binary systems. We discussed 
the code in detail in Paper I, but in short, it utilizes 
the Eclipsing Light Curve (ELC) code, which was writ- 
ten to compute light curves of eclipsing binary systems 
(Orosz & Hauschildt 2000). The ELC code represents 
the surfaces of two stars as a grid of individual lumi- 
nosity points, and calculates the resulting light curve 
given the provided systemic parameters. ELC includes 
the dominant physical effects that shape a binary's light 
curve, such as non-spherical geometry due to rotation, 
gravity darkening, limb darkening, mutual heating, re- 
flection effects, and the inclusion of hot or cool spots on 
the stellar surface. For the work in this paper we have 
simply turned off one of the stars, thus allowing us to 
probe the astrometric effects of a single star. To com- 
pute intensity, ELC can either use a blackbody formula 



or interpolate from a large grid of NextCen model atmo- 
spheres (Hauschildt et al. 1999). For all the simulations 
in this paper, we have used the model atmosphere op- 
tion, and will note now, and discuss more in detail later, 
that the calculation of limb-darkening is automatically 
included in NextCen model atmospheres. These arti- 
ficially derived limb-darkening coefficients have recently 
been shown to be in error by as much as ~1 0-20% i n com- 
parison to observationally derived values ()Claretl [2008'). 
and thus their uncertainties must be included, although 
for this work, due to symmetry, we find the introduced 
error is negligible. For all our simulations, we model the 
U, B, V, R, I, J, H, and K bands for completeness and 
comparison to future studies, though we note that SIM 
Lite will not be able to observe in the U, J, H, or K 
bandpasses. 

3. INCLINATION AND ROTATION 

The astrophysical phenomenon of gravity darkening, 
also sometimes referred to as gravity brightening, is the 
driving force behind the ability to determine the inclina- 
tion of a single star using multi-wavelength astrometry. 
A rotating star is geometrically distorted into an oblate 
spheroid, such that its equatorial radius is greater than 
its polar radius, and thus the poles have a higher surface 
gravity, and the equator a lower surface gravity, than a 
non-rotating star with the same mass and average radius. 
This increased surface gravity, g, at the poles results in a 
higher effective temperature, T^s, and thus luminosity; 
decreased g at the equator results in a lower Teg and 
luminosity. This temperature and luminosity differen- 
tial causes the star's center of light, or photocenter, to 
be shifted towards the visible pole, away from the star's 
gravitational center of mass. Since the inclination de- 
termines how much of the pole is visible, the amount of 
displacement between the photocenter and the center of 
mass is directly related to the inclination. Furthermore, 
since the luminosity difference effectively results from a 
ratio of blackbody luminosities of differing temperatures, 
the effect is wavelength dependent, with shorter wave- 
lengths shifted more than longer wavelengths. Thus, the 
amount of displacement between the measured photocen- 
ter in two or more wavelengths is directly related to the 
inclination. See Figure [T] for an illustration of the effect. 

An additional complicating factor i s the exact depen - 
dence of temperature on local gravity. Ivon Zeipell (IT92I 
was the first to derive the quantitative relationship be- 
tween them, showing that T^g oc g^'^ , where /3i is referred 
to as the gravity darkening exponent. The value of /3i 
has been a subjec t of mu c h stud y and debate; for a com- 
plete review, see iClareti (|200 0f). who presents both an 
excellent discussion of past studies, as well as new, de- 
tailed computations of Pi using modern models of stellar 
atmospheres and internal structure that encompass stars 
from 0.08 to 40 Mq. Since the value of (3i affects the 
temperature differential between equator and pole, the 
multi- wavelength displacement will also be dependent on 
the value of (3i . The total amplitude of the effect will be 
scaled by the angular size of the star, which depends on 
both its effective radius and distance. Thus, in total, 

^ REFLU X can be run via a web interface from 
http : //astronomy . nmsu . edu/j Icough/ref lux . html Additional 
details as to how to set-up a model arc presented there. 
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Figure 1. Illustration of the inclination effect using brightness maps of Capella Ab, which has an inclination of 42.788° (To rres et al.l 
12003) . Left: Capella Ab artificially spun- up to near its break-up speed, to accentuate the gravity darkening effect for ease of viewing. As 
can be seen, the photocenter of the system is dramatically shifted away from the center of mass, towards the visible pole, which is brighter 
than the rest of the star due to gravity darkening. Furthermore, since the pole is physically hotter, the U-band photocenter is shifted 
more than the K-band photocenter, and in the direction of the projected rotation axis, or y-axis. Right: Capella Ab at its actual rotation 
period. As can be seen, the actual effect is small compared to the angular size of the star on the sky, but still large compared to the 1 fias 
benchmark of SIM Lite. The presence of gravity darkening is clearly visible, which causes a decrement in flux towards the limb of the star. 
Note that the broad wavelength covereage of SIM Lite will only cover the B, V, R, and I bandpasses. 



the components of this inchnation effect are the effective 
stellar radius, distance, effective temperature, rotation 
rate, and inclination of the star. In principle, one 
is able to determine the effective stellar radius, effective 
temperature, rotation rate, and distance of a target star 
using ground-based spectroscopy and space-based paral- 
lax measurements, including from SIM Lite. Thus, when 
modeling the multi-wavelength displacement of the stel- 
lar photocenter, the only two components that need to 
be solved for are the inclination and /3i, with /3i already 
having some constraints from theory. 

A good trio of stars for modeling and testing this in- 
clination effect are the components of the binary sys- 
tem Capella , (Aa and Ab), and the single star Vega. 
iTorres et al.l (|2009l ) has very recently published an ex- 
tremely detailed analysis of both the binary orbit of 
Capella and the physical and evolutionary states of the 
individual components, providing both new observations, 
as well as drawing from the pre vious obs ervatio ns and 
analys es of iHummel et al.l ()1994l ) and iStr assmeie r et all 
(IMI . Ve ga, in addition to being one of the most 
well-studied stars in the sky, has recently been discov- 
ered to be a very rap i d rotator seen near l y pole-on 
(lAufdenberg et al.H2006l: iPeterson et"aLl [20061: lllill et al.l 
|2010[) . In total, these three stars represent both slow 
and rapid rotators for giant and main-sequence stars at 
a range of temperatures, as Capella Aa is a slow-rotating 
K-type giant, Capella Ab is a fast- rotating G-type giant, 
and Vega is a very fast-rotating A-type main-sequence 
star. With many ground-based interferometric obser- 
vations to compare with, and being bright and nearby, 
these stars also present excellent targets for SIM Lite. 

We use the reflux code to generate models of the 
astrometric displacement from U -band to 77-band, with 
respect to the i^-band photocenter, for inclinations from 



Table 1 

Parameters for Capella Aa 



Parameter 


Value'' 


Distance (pc) 


12.9 


Rotation Period (Days) 


106.0 


Mass (Mq) 


2.70 


Radius (R©) 


12.2 


Effective Temperature (K) 


4940 


/3i 


0.43 



=■ Values from ITorres et al.l 120091') 



to 90°, for each star, as shown in Figures [51 ^ and Bl 
We use systemic parameters giv en by [T orres ct al. (2009^ 
for Capel la Aa and Ab, a nd by lAufdenberg et al.i ((2006r ) 
and Peter son et all (j2006[) for Vega, listed in Tables [Tl[2l 
and [3| respectively. We employ the model atmospheres 
incorporated into the ELC code, as well as aut omati- 
cally chosen values for /3i based on Figure 1 of IClaretl 
(|2000D . Additionally, in each figure we show a dashed 
line to indicate the effect of decreasing the gravity dark- 
ening coefficient by 10% to simulate the uncertainty of 
the models (Claret 2000) and explore the correlation with 
other parameters. 

As can be seen from these models, we find that the 
effect is quite large for a Capella Ab-like or Vega-like 
fast rotator, but only marginally detectable for a slower- 
rotating system like Capella Aa. This also implies that 
this effect would not be detectable for a slow-rotating, 
main-sequence star like our Sun. Our modeling confirms 
this, showing a total U-K amplitude of ^ 0.1 /las for 
a 1.0 Mq, 1.0 Rq star with a rotation period of 30.0 
days at 10.0 parsecs. These conclusions on detectability 
are made with the assumption that, for bright stars like 
these, SIM Lite can achieve its microarcsecond bench- 
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Figure 2. The astrometric displacement of each bandpass with 
respect to /f-band versus inclination for Capella Aa. Dashed lines 
are a model with fi\ decreased by 10%. Due to the slow rotation 
rate of Capella Aa, the effect is limited to a maximum of ~2.0 
microarcseconds between U and K-band, and only a maximum of 
'^0.7 /jas between B and I-band, where SIM Lite will operate. This 
puts the detection of this effect for Capella Aa at the very edge of 
SIM Lite's capability. 




Inclination (Degrees) 

Figure 3. The astrometric displacement of each bandpass with 
respect to K-band versus inclination for Capella Ab. Dashed lines 
are a model with fi\ decreased by 10%. Due to the fast rotation 
rate of a Capella Ab like star, the effect is moderate with tens of 
microarcseconds of displacement, and thus these types of stars are 
excellent targets for SIM Lite. The actual inclination of Capella Ab 
is 42.788° (Torr es et al.,,200 9'). and thus Capella Ab itself should 
show a large shift of the photocentcr with wavelength. Note that 
the broad wavelength covereage of SIM Lite will only cover the B, 
V, R, and I bandpasses. 



Table 2 

Parameters for Capella Ab 



Parameter 


Value'' 


Distance (pc) 


12.9 


Rotation Period (Days) 


8.64 


Mass (Mq) 


2.56 


Radius (Rq) 


9.2 


Effective Temperature (K) 


5700 


/3i 


0.39 


Values fromlTorres et al.l 1200911 
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Figure 4. The astrometric displacement of each bandpass with 
respect to /f-band versus inclination for Vega. Dashed lines are a 
model with /3i decreased by 10%. Note that due to the very fast 
rotation of Vega, along with a high value of /3i , the effect can be 
quite large, at a couple hundred micro-arcseconds. For a Vega-like 
star, SIM Lite observations would yield very accurate values for /3i 
and the inclination. For Vega itself , which is know n to be nearly 
pole on, with an inclination of 5.7° IIHill et al.ll20Tot) . there should 
be a B-band minus I-band displacement of 6.0 /xas, still detectable 
by SIM Lite. Note that the broad wavelength covereage of SIM 
Lite will only cover the B, V, R, and I bandpasses. 



Table 3 

Parameters for Vega 



Parameter 


Value" 


Distance (pc) 


7.76 


Rotation Period (Days) 


0.521 


Mass (M0) 


2.11 


Radius (R©) 


2.5 


Effective Temperature (K) 


9602 


/3i 


1.02 



" Va lues from j Aufdenbcrg et al] 
pOM ) and lPeterson et al.. (,2006 ) 



mark. We show this is possible in narrow angle (NA) 
mode by employing the SIM Dif ferential Astrorn etry Per- 
formance Estimator (DAPE) (jPlummeii |2009[) . For a 
target star with magnitude V=5, and a single compar- 
ison star with V=10 located within a degree of it on 
the sky, by integrating 15 seconds on the target, and 
30 seconds on the reference, for 10 visits at 5 chop cy- 
cles each, a final precision of ±1.01 /xas is achieved in 
only 1.04 hours of total mission time. For a fainter tar- 
get with T^=10, this precision is only reduced to ±1.32 
/ias in the same amount of mission time. In utilizing 
NA mode, one must be careful in choosing the refer- 
ence star(s), to ensure that they are not stars with a 
substantial wavelength dependant centroid. Given the 
only constraints on reference stars are that they need to 
have V > 10 and are within one degree on the sky, one 
could easily choose a slow-rotating, main-sequence star, 
determined as such via ground-based observations, as a 
wavelength-independent astrometric reference star. We 
also note that wide angle SIM Lite measurements, with 
a precision of ~5 /.tas, may not detect the wavelength 
dependent photoceter of a system like Capella, but will 
have no difficulty detecting it in stars like Capella Ab or 
Vega. 
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The effect of decreasing the gravity darkening exponent 
is to decrease the total amphtude of the effect in each 
wavelength, with shorter wavelengths affected more than 
longer wavelengths. Thus, the choice of gravity darken- 
ing exponent is intimately tied to the derived inclination. 
If one were to model observed data with a gravity dark- 
ening exponent that was ~10% different than the true 
value, they would derive an inclination that would also 
be ~10% different from the true inclination. However, 
the two combinations of inclination and gravity darken- 
ing exponent do not produce identical results, and can 
be distinguished with a sufficient precision at a number 
of wavelengths. For example, if one were to adopt the 
nominal value for Pi and derive an inclination of 40 de- 
grees for a Vega-like star, then adopt a f3i value that was 
10% lower, one would derive an inclination of 43 degrees, 
a 7.5% change. In this case though, with the lower /3i 
value, the measured photocenter in the U, B, V, R, I, 
J, and H bandpasses, with respect to the K-band photo- 
center, would differ from the nominal /3i model by ~0.5, 
-1.0, -2.0, -2.0, -1.6, -1.0, and 0.2 /xas respectively. Note 
that for B, V, R, and I, where SIM Lite can observe, 
these discrepancies, on the order of '-^1.0 /xas, should be 
large enough to be distinguished in NA mode. Thus, a 
unique solution exists for the values of i and /3i if the 
photocenter is measured in three or more wavelengths. 
(The photocenter of one wavelength is used as a base 
measurement that the photocenters of other wavelengths 
are measured with respect to, as we have chosen K-band 
as the base measurement in our models. With the pho- 
tocenter measured in three or more wavelengths total, 
there are two or more photocenter difference measure- 
ments, with two unknown variables for which to solve.) 
Another complication is the possibility of having equally 
good fitting high and low solutions for i. For example, if 
one observed and determined a best-fit inclination of 70 
degrees for a Vega-like star, one could obtain a reason- 
ably good fit as well at 46 degrees, (see FigS]). However, 
just as in the case of the uncertainity in the value of /3i, 
discernible discrepancies would exist. In this case, the 
discrepencies in the measured photocenter in the U, B, 
V, R, I, J, and H bandpasses, with respect to the K-band 
photocenter, would be ~0.1, -9.0, -2.0, 1.5, 6.0, 1.0, and 
0.2 /xas respectively. Just as in the case of the uncertain- 
ity in the value of this discrepancy between equally 
good fitting high and low inclination solutions can be re- 
solved if one has three or more wavelengths obtained in 
NA mode. 

As mentioned in [JTJ we note that the limb- darkening 
function, which was automatically chosen by the ELC 
code as incorporated into the model atmospheres, can 
differ from actual observed values by ^-^10% (Claret 
12008*) . We have tested how changing the limb-darkening 
coefficients by 10% affects the resulting astrometric dis- 
placements, and find that the result is less than 0.5% for 
all wavelengths, and thus is negligible in the modeling. 
The reason is that limb-darkening is symmetric, and thus 
while increased limb-darkening damps the visible pole, it 
also damps the rest of the star, and thus the relative 
brightness between regions is maintained. 

Additionally, this inclination technique yields the ori- 
entation of the projected stellar rotation axis on the sky, 
which is parallel to the wavelength dispersion direction. 
When coupled with the derived inclination, this tech- 



nique thus yields the full 3-dimensional orientation of 
the rotation axis. This could be a powerful tool in deter- 
mining the overall alignment of stellar axes in the local 
neighborhood and in nearby clusters. 

4. STAR SPOTS 

Another area of astrophysical interest to which multi- 
wavelength astrometric measurements from SIM Lite can 
contribute is the study of star spots. As the cause of 
star spots are intense magnetic fields at the photosphere, 
they are typically found in stars with convective en- 
velopes, especially rapidly rotating stars. Thus, both 
low-mass, main-sequence K and M dwarfs, as well as 
rapidly rotating giant and sub-giant stars, are known 
to host large spots on their surface. The study of 
the distribution, relative temperature, and size of these 
spots would greatly contribute to the study of mag- 
netic field generation in stellar envelopes. A starspot 
that rotates in and out of view will cause a shift 
of the photocenter for a single star, which has been 
a su bject of much recent discu s sion in the literature 
(e.g. Hatzcs 2002; Unwin 2005; Eriksson & Lindcgren| 
2007; Catanzaritc et al. 2008; Makarov ct al. 200^ 
Lanza et al.ll2008t ). especially in light of its potential to 
mimic, or introduce noise when characterizing, an ex- 
trasolar planet. However, there has been no mention in 
the literature of the multi- wavelength astrometric signa- 
ture of stellar spots, where, just as in the case of the 
gravity darkening inclination effect, we are looking at es- 
sentially two blackbodies with varying temperatures, and 
thus shorter wavelengths will be more affected by a spot 
than longer wavelengths. 

To characterize the multi-wavelength astrometric sig- 
nature of stellar spots, we model two spotty systems, 
again using the reflux code. We model Capella Ab, 
which shows evidence of la rge spots and is susp ected of 
being a RS CVn variable ([Hummel et al.|[l994[ ). and a 
typical main-sequence K dwarf. For Capella Ab, we use 
the parameters listed in Tabled along with the star's de- 
termined inchnation of 42.788° (Torres et al. 2009), and 
add a cool spot that has a temperature that is 60% of the 
average surface temperature, located at the equator, at a 
longitude such that it is seen directly at phase 270°, and 
having an angular size of 10°, (where 90° would cover 
exactly one half of the star). For the K dwarf system, 
we use the physical parameters listed in Table |4l simu- 
lating a typical K Dwarf at 10 parsecs, and add a cool 
spot with the same parameters as we do for Capella Ab. 
Additionally, to investigate the effects of cool versus hot 
spots or flares, we also run a model with a hot spot by 
changing the spot temperature to be 40% greater than 
the average surface temperature. We present our models 
in Figures [5l El and [71 

As can be seen for CapellaAb, the gravity darkening 
inclination effect presented in ^ dominates the spread 
of colors in the y-direction, the direction parallel to the 
stars' projected rotation axis. However, the amplitude 
of the spot motion is quite large, with a total ampli- 
tude of '-~^40 /ias in all bandpasses, which would be easily 
detectable by SIM Lite. For the K dwarf with a cool 
spot, we see a much smaller, but still detectable shift of 
amplitude '-^5-8 /zas, depending on the wavelength. In 
the case of a hot spot or flare, we see a much larger 
displacement, on the order of ~10-200 //as, depending 
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Figure 5. A simulated cool Spot on CapellaAb. The spot is located on the equator, with a longitude such that it is seen directly at phase 
270°. The strong presence of the gravity darkening effect, discussed in [[3] dominates the wavelength spread in the y direction. Note that 
the broad wavelength covereage of SIM Lite will only cover the B, V, R, and I bandpasses. 
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X (nas) Phase (degrees) 

Figure 6. A simulated cool spot on a nearby K dwarf star with an inclination of 60°, whose parameters are given in Table |4] The spot 
is located on the equator, with a longitude such that it is seen directly at phase 270° . Note that the broad wavelength covereage of SIM 
Lite will only cover the B, V, R, and I bandpasses. 
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X (lias) Phase (degrees) 

Figure 7. A simulated hot spot or flare on a nearby K dwarf star with an inclination of 60°, whose parameters are given in Tabled The 
spot is located on the equator, with a longitude such that it is seen directly at phase 270°. Note that the broad wavelength covereage of 
SIM Lite will only cover the B, V, R, and I bandpasses. 
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Table 4 

Parameters for the K Dwarf System 



Parameter 


Value 


Distance (pc) 


10.0 


Inclination (°) 


60.0 


Period (Days) 


20 


Mass (M0) 


0.6 


Radius (R©) 


0.6 


Effective Temperature (K) 


4500 


Latitude of Spot (°) 


90 


Longitude of Spot (°) 


270 


Angular Size of Spot (°) 


10 


Cool Spot Temperature Factor 


0.6 


Hot Spot Temperature Factor 


1.4 



on the wavelength, which would be easily detectable by 
SIM and provide extremely precise values in deriving the 
spot parameters. 

In general, the temperature of the spot, in relation to 
the mean stellar surface temperature, is related to the 
spread in observed wavelengths, with a larger spread in- 
dicating a larger temperature difference. The duration of 
the astrometric displacement in phase, coupled with the 
overall amplitude of the astrometric displacement, yields 
the size of the spot, as larger spots will cause larger dis- 
placements and be visible for a larger amount of rota- 
tional phase. The latitude of the spot can also affect the 
total duration. Finally, the amplitude of the astrometric 
displacement in the x versus the y direction is dependent 
on both the latitude of spot as well as the inclination of 
the star. Thus, when modeled together, one is able to 
recover these parameters. This work can also be com- 
bined with our work in Paper I to derive the location of 
spots in binary systems, as the astrometric signature of 
the spot is simply added to the astrometric signature of 
the binary system. 

The astrometric motion induced upon a parent star by 
a host planet does not have a wavelength dependence. 
Spots however, as we have shown via our modeling, have 
a clear wavelength dependence. Thus, if one has a can- 
didate planetary signal from astrometry, but it shows a 
wavelength-dependant motion, it must then be a false 
positive introduced from star spots at the rotation pe- 
riod of the star, (assuming that the planet's emitted flux 
is negligible compared to the star.) Furthermore, when 
SIM is launched, there will likely be many cases where 
a marginally detectable signal due to a planetary com- 
panion is found at a very different period than the ro- 
tation period of the star. However, starspots will still 
introduce extra astrometric jitter which will degrade the 
signal from the planetary companion. Multi-wavelength 
astrometric data can be used to model and remove the 
spots, which will have a wavelength dependence, and 
thus strengthen the planetary signal, which will not have 
a wavelength dependence. 

5. DISCUSSION AND CONCLUSION 

We have presented detailed models of the multi- 
wavelength astrometric displacement that SIM Lite will 

observe due to gravity darkening and stellar spots using 
the REFLUX code. We find that SIM Lite observations, 
especially when combined with other techniques, will be 
able to determine the absolute inclination, gravity dark- 
ening exponent, and 3-dimensional orientation of the ro- 



tational axis for fast and slow rotating giant stars, and 
fast-rotating main-sequence stars. This technique will 
be especially useful in probing binary star and exoplanet 
formation and evolution, as well as the physics of star 
forming regions. Direct observational determination of 
the gravity darkening exponent has direct applications in 
both stellar and exoplanet astrophysics. This technique 
is also relatively inexpensive in terms of SIM Lite observ- 
ing time, as one need only to observe a given star once, 
as opposed to binary stars and planets, which require 
constant monitoring over an entire orbit. It should be 
noted that this effect should be taken into account when 
constructing the SIM Lite astrometric reference frame, 
such that fast-rotating giants should be cxcJudcid so as 
not to produce a wavelength-dependent astrometric ref- 
erence fame. 

We also have presented models of star spots on sin- 
gle stars, and find that SIM Lite should be able to dis- 
cern their location, temperature, and size. Combined 
with other techniques, this will provide great insight into 
stellar differential rotation, magnetic cycles and underly- 
ing dynamos, and magnetic interaction in close binaries. 
From this modeling, it should especially be noted that 
multi- wavelength astrometry is a key tool in the hunt for 
extrasolar planets, either by ruling out false signals cre- 
ated by spots, or simply removing extra astrometric jitter 
introduced by spots. Thus, it remains critical that SIM 
Lite maintains a multi- wavelength astrometric capability 
in its final design. 
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Study (PI: D. Gelino) from the National Aeronautics and 
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Propulsion Laboratory, California Institute of Technol- 
ogy. J.L.C. acknowledges additional support from a New 
Mexico Space Grant Consortium Fellowship. We thank 
the referee for comments which greatly helped to improve 
this manuscript, particularly in making the presentation 
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